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ABSTRACT 

It is expected on both evolutionary and empirical grounds that many merging neutron star (NS) 
binaries are composed of a highly magnetized NS in orbit with a relatively low magnetic field NS. I 
study the magnetic interactions of these binaries using the framework of a unipolar inductor model. 
The e.m.f. generated across the non-magnetic NS as it moves through the magnetosphere sets up a 
circuit connecting the two stars. The exact features of this circuit depend on the uncertain resistance 
in the space between the stars 7\L S p ace . Nevertheless, I show that there are interesting observational 
and/or dynamical effects irrespective of its exact value. When lZ apacc is large, electric dissipation 
as great as ~ 10 46 erg s _1 (for magnetar-strength fields) occurs in the magnetosphere, which would 
exhibit itself as a hard X-ray precursor in the seconds leading up to merger. With less certainty, 
there may also be an associated radio transient, but this would be observed well past merger (~ hrs) 
because of interstellar dispersion. When lZ spacc is small, electric dissipation largely occurs in the 
surface layers of the magnetic NS. This can reach ~ 10 49 erg s _1 during the final ~ 1 s before merger, 
similar to the energetics and timescales of short gamma-ray bursts. In addition, for dipole fields 
greater than ss 10 12 G and a small TZ space , magnetic torques spin up the magnetized NS. This drains 
angular momentum from the binary and accelerates the inspiral. A faster coalescence results in less 
orbits occurring before merger, which would impact matched-filtering gravitational-wave searches by 
ground-based laser interferometers and could create difficulties for studying alternative theories of 
gravity with compact inspirals. 

Subject headings: binaries: close — gamma rays: bursts — gravitational waves — stars: magnetic 
fields — stars: neutron 
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1. INTRODUCTION 

The inspiral and coalescence of a neutron star (NS) 
with another NS or a black hole is generally consid- 
ered the most promising target for ground-based grav- 
itational wave det ectors. It is anticipated t hat the inter- 
ferom eters LI GO ([Abramovici et al.lll992t lAbbott et al.l 
120091) . Virgo (iCaron et alj Il999t lAcernese et all 12009). 
and KAGRA (|Kuroda et al.l 120101 ) will reach sufficient 
sensitivity to observe these events within the next few 
year s once they are up graded to "advanced" sensitiv- 
ity (jAbadie et all 120101 ). Since tidal interactions can- 
not enforce synchroniza tion in a coalescing NS binary 
(Bilds ten fc Cutlerlll992D . the associated waveforms are 
expected to be relatively simple and can be predicted 
with high accuracy for the inspiral phase. This makes 
them ideal for using matched-filteri ng methods for ex- 
tract ing their signal from the noise ([Cutler fc Flanaganl 
1994). In addition, these merger events are of inter- 
est for their potential electromagnetic signatures. They 
are the favored progenitors of short gamma-ray bursts 
CGRBs. IBlinnikov et al.l [l98l IPaczvnskil [19861 fl99Tt 
lEichler et al.l 119891 ) and have been theorized to produce 
~ 1 day timescale opti cal transients from the eje ction of 
neutron-rich material (|Metzger et al.ll2008l l2010h . 

Strong magnetic fields are undoubtedly present in a 
large fraction of these merging systems. In particu- 
lar, it is expected that in many binaries one NS has 
a high magnetic field (~ 10 12 — 10 14 G) and the other 
has a lower-strength field (~ 10 8 — 10 9 G). Such 
a situation is anticipated on theoretical grounds from 



the binary evolution of two massive m ain-sequence star 
(|Bhattacharva fc van den Heuve]lH99l . The more mas- 
sive star experiences a supernova first and forms a NS. 
When the less massive star reaches the end of its main 
sequence, it becomes a red supergiant before going super- 
nova and producing a NS of its own. This transfers an- 
gular momentum a nd mass to the first NS and decreases 
its magnetic field (iBisnovatvi-Kogan fc Komberel H974t 
iShibazaki et all 119891) . The result is an older NS with 
a magnetic field similar to a recycled millisecond pulsar 
in a binary with a younger NS that has a magnetic field 
more typical of a normal radio pulsar. Such a picture has 
been confirmed in at least one case fo r the double pulsa r 
PSR J0737-3039 (IBurgav et al.l 120031 ILvne et alll2Ciol . 
in which the fast and slow spinning NSs have dipole fields 
of 6.3 x 10 9 G and 1.2 x 10 12 G, respectively. In addition, 
since ~ 10% of pulsars ar e born with magnetar strengt h 
fields (~ 10 14 - 10 15 G; IKulkarni fc Thompson! [l998l) . 
there is the potential for some binary NSs to have ex- 
tremely strong fields present (although this may depend 
on what evolutionary scenarios are needed to produce 
such high magnetic fields). 

The effect of magnetic fields on NS coalescence has 
been considered in a number previous theoretical stud- 
ies. Numerical simulations have investigated how the 
fluid dynamics of the NSs are i mpacted by magnetic 
fields, especiall y post merger (see lAnderson et al.l 120081 : 
iLiu et all 120081: iGiacomazzo et ali 120091 . 12011L and ref- 
erences therein). Other work has focused on the in- 
spiral phase to understand the interaction of a non- 
magnetized NS moving through another NS's mag- 
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netosphere (iLipunoy fe Panchenkol 119961 : iVietril 119961 : 
Hansen fc Lvutikovl l200lT T My work here has a focus 
similar to these latter studies, but with greater empha- 
sis on the time-dependent binary evolution. The asyn- 
chronicity of stars during inspiral provides free energy 
that can be tapped to drive currents and electrical dissi- 
pation. The strength and location of this dissipation is 
strongly influenced by the (uncertain) resistivity of the 
space between the merging NSs, 7^ spacc . Even with this 
large uncertainty, I show that interesting observational 
and/or dynamical effects occur no matter what value of 
7?- S paco is chosen. 

In 32 I describe the main features of the so-called 
"unipolar inductor" model applied to binary NSs. In [J3l 
I investigate the interplay between magnetic torques and 
the loss of angular momentum due to gravitational wave 
emission. I provide analytic estimates for the resulting 
level of synchronization and the electrical power dissipa- 
tion as a function of orbital frequency. In SjH I integrate 
the time evolution of the coalescing binary numerically 
to solve for the effects of magnetic fields. I especially fo- 
cus on how the dissipation in the circuit connecting the 
two NSs is impacted by 72. spacc . In fJSJ I conclude with 
a summary of my main results, provide a discussion of 
the potential electromagnetic signatures of these mag- 
netic interactions, and highlight many of the outstand- 
ing problems that remain in understanding the impact 
of magnetic fields on inspiralling NSs. 

2. THE UNIPOLAR INDUCTOR MODEL 

Consider a NS-NS binary merging from gravitational 
wave emission. One NS has a mass M\ and radius R\ and 
the other has mass M 2 and radius R 2 . When the binary 
has an separation a, corresponding to a frequency w, the 
orbital energy is 



£orb — 



GM X M 2 



= (Geo) 2 / 3 



M X M 2 
MV3 : 



(1) 



where M = M\ + M 2 . (For this paper I ignore post- 
Newtonian corrections since they are not crucial to any 
of the main results.) At a gravitational wave frequency 
/ = cu/ir = 10 2 Hz, which is near the optimum region 
for detection by ground-based laser interferometers, the 
energy is £" rb ~ 3 x 10 52 erg. In contrast, the rotational 
energy of Mi when it has spin Oi is 



E IQ t — 2^i^i' 



(2) 



where I\ = 0.35MjR 2 is the app roximate moment of in- 
ertia (|Lattimer fe P rakash 2001). In the extreme limit of 
being synchronized, Oi « u> and E YOt m 6 x 10 49 erg. In 
other words, if merely 0.2% of the orbital energy is put 
into a NS's spin, this is enough for the NS to be synchro- 
nized with the orbit. Furthermore, typical energies for 
the prompt emission from short GRBs (once corrected for 
beaming) are - 10 48 - 10 50 erg (|Nakarll2007lh So again, 
if only a small fraction of the orbital energy is tapped 
(and then emitted on a sufficiently short timescale), it 
may result in an interesting observable signal. In the fol- 
lowing work I argue that magnetic effects may provide 
the coupling to tap this energy. 

For the majority of this paper I use the general frame- 
work of a unipolar inductor model to assess the interac- 
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Fig. 1. — Schematic diagram summarizing the main features of 
the unipolar inductor model applied to binary NSs. The mag- 
netic NS has a dipole field, which threads the non-magnetic NS. 
The orbital frequency is ui, and the magnetic NS has a spin fre- 
quency f2i . Therefore, the non-magnetic NS has a relative velocity 
v with respect to the magnetosphere, which points into the page 
when w > Ct\. In the frame of the non-magnetic NS, there is an 
induced electric field E = (vxB)/c (the magnetic field roughly 
points downward at this location), which drives a current loop I 
that is traced out with arrows. For the majority of the loop, the 
current is forced to flow along the magnetic field lines, with the 
exception of the surfaces of the two NSs. The interaction of the 
current over the non-magnetic NS with the dipole field drives a 
magnetic toque = 2i?2[rX (XxB)/c\. This points in the same 
direction as f2i when u > tli, spinning up the magnetic NS. This 
diagram only focuses on the upper half of the magnetic interac- 
tions. A similar process is repeated in the lower hemisphere which 
adds to the torque in the same direction. 

tion of a non-magnetic NS with a magnetic NS. The main 
ideas of this model have been presented and explored in 
a variety of astrophysical setti ngs in the past, from the 
inter action of Jupiter and Io (|Goldreich fc Lynden-Belil 
1969) , to ultracompact whi te dwarf binaries (|Wu et al.l 
120021; IPairOsso et al.ll2006fl. to sho rt orbital period ex- 
trasolar Earths ( Laine fc Linl I2012J). to the interaction 
of a magnetic NS and black hole (|Mc Willia ms fc Levinl 
120111 ). Below I discuss this model in the context of 
merging NSs. The main ideas are summarized schemat- 
ically in Figures Q] and [21 and I will be referring to these 
throughout the discussion. 

2.1. Skin Depth and Voltage Estimates 

One NS (labeled as "1") has a magnetic moment 
/i = BR 3 , where B is its dipole magnetic field strength. I 
assume that the magnetic moment, the spin axis of each 
NS, and the spin axis of the orbit are all aligned, as shown 
in Figure [T] The other NS (labeled as "2" ) has a neg- 
ligible magnetic field. The ability of the magnetic field 
to diffuse into the non-magnetic NS is determined by its 
cond uctivity, which is roughly given by (|Spitzer fc Harm! 
1953) 



K = 7 



2(2fc B r e /7r) 3 / 2 

m!' 2 Ze 2 In A 



(3) 



where fce is Boltzmann's constant, T e is the electron tem- 
perature, m e is the electron mass, e is the electron charge, 
Z is the charge per ion, and In A is the Coulomb loga- 
rithm. The factor 7 depends on the charge per ion Z, 
and can vary between 0.6 (for Z = 1) and 1 (in the limit 
Z — > 00). The characteristic conductivity is therefore 
if a 1.4 x 10 16 s" 1 , where I take T e w 10 6 K, Z = 1, 
7 « 0.6, and In Aw 10. 

Since the binary is inspiralling due to gravitational 
wave emission, there is only a short amount of time for 
magnetic diffusion to occur. The rate of angular momen- 
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turn l oss due to gravitational waves is (jLandau fc Lifshit zl 
1975) 



32 G 3 M X M 2 M 



Jo 



(4) 



where J or b = (Ga/M) 1 / 2 MiM2- The inspiral timescale 
is 



J, 



orb 



M 1 / 3 



31 J, 



96 G 5 / 3 MiM 2 w 8 / 3 



140 



10 2 s- 



-8/3 



(5) 



where I have assumed that each NS has a mass of 1.3M©, 
and the factor of 1/3 is because I use the timescale 
associated with how the orbital frequency is changing 
(J or b oc w -1 / 3 ). As described in lJacksonl ([1975J), the skin 
depth for this timescale is roughly 



2nK 



1/2 



which as a ratio to the non-magnetic NS radius is 



R 2 



1.2 x lO -3 ^- 



r 6 1/2 ( 



10 2 



-4/3 



(6) 



(7) 



where Ki§ = K/1Q 16 s 1 and I have used a radius of 
R 2 = 12 km. 

Since this skin depth is small in comparison to the 
radius, I assume that the non-magnetic NS is roughly 
a perfect conductor . As d escribed in the Appendix of 
lHansen fc Lyutikovl (|2001[ ). in this limit the magnetic 
field must be completely excluded from the non-magnetic 
NS. This implies that surface currents are induced which 
produce a magnetic dipole with an opposite orientation 
with respect to the downward pointing dipole field at 
the location of the non- magnetic NS. If the non- magnetic 
NS is also spinning, one must take into account this to- 
tal magnetic field configuration to correctly understand 
the interaction. Since the effect of the non- magnetic NS's 
spin is subdominant to the orbit, and it has alre ady been 
well-summarized by lHansen fc Lyutikovl (|2001[ ) , I ignore 
it here. 

Instead, for this study I focus on the orbital motion of 
the non-magnetic NS. This induces a surface charge den- 
sity with a dipole structure pointed toward the magnetic 
NS 



S = B (e r xv)/4nc, 



(8) 



where e r is a radial unit vector with respect to the cen- 
ter of the non-magnetic NS, and v is the non- magnetic 
NS's velocity relative to the magnetosphere. Due to this 
charge density distribution, there is a net electric field 
from one side of the non-magnetic NS to the other of 
E = (vxB)/c, and an associated potential difference 
$ ss 2R,2\E\. The e.m.f. when the NSs are at a separa- 
tion a is therefore 



4> 



2^R 2 



(lu - fix) 



7/3 



01 



c(GM) 2 / 3 u ' 



(9) 
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Fig. 2. — Diagram of the circuit formed between the two NSs. 
The non-magnetic NS has an induced e.m.f. <E>, which drives a 
current I over the circuit. Each leg of the circuit has an associated 
resistivity, which are denoted as 7^i , 7^2 1 and 7?. S pacc ■ The relative 
size of each of these resistances determine the location and amount 
of dissipation. 

the magnetic NS. The numerical value for the voltage is 

x 1013 m Go 2 !- 1 ) 773 (SO statvoit ' ( io ) 

where ^31 = ^/10 31 G cm 3 . Note that this strength 
magnetic moment corresponds to a magnetic field of 
w 6 x 10 12 G, which highlights the fact that a surpris- 
ingly large potential is possible even for a magnetic field 
much lower than that typically associated with magne- 
tars. 

2.2. Locations and Rates of Dissipation 

The induced electric field across the non-magnetic NS 
drives a current over a circuit that connects the two stars, 
as shown in Figure [2] For a given circuit, the sources of 
resistance are: the magnetic NS, the non-magnetic NS, 
and the two regions of space connecting the stars, which 
are denoted as TZ\, H2, and 7?. spacc , respectively. Since 
the resistors are simply connected in series, the current 
is I = <&/{1Zi + IZ2 + 27vL spaC c)- Through each resistor 
there is an associated rate of dissipation. For example 
through the magnetic NS, 



Ei = 2$ 2 ^x/(^i +TZ 2 + 2TZ space ) 



(11) 

and likewise for each of the other lengths of the cir- 
cuit (note a factor of two is included because the cir- 
cuit is repeated in both hemispheres). From this one can 
see that the relative sizes of the resistances will have a 
big influence over where the dissipation is greatest and 
by what amount. For example, if IZi 3> IZ2 j 72-spacc , 
then Ei « 2$ 2 /^i, but if ft spacc > n 1 ,n 2l then 
Ei^^TZ 1 /2TZ% acc ^2^/n 1 . 

For a characteristic vertical length L and area A, 
the resistance is related to the conductivity K by 
1Z = L/AK. Now consider the magnetic field lines that 
simultaneously thread both NSs. For the non- magnetic 
NS, current flows over its entire surface. But when the 
connecting field lines are traced back to the magnetic 
NS, they only intersect a small fraction of its surface 
area (see Figure [T]). Due to this difference in intersected 
areas, I approximate that IZi 3> TL% and thus ignore IZ2 
in the rest of my study. Using the geometry of the binary 
and dipole field, IWu et al. I ([2001 estimate 



Tli 



J 
R2K 



Ad J 



a 



3/2 



(12) 



where ai = 2(oj — Oi) is the tidal forcing frequency on where J is a geometric factor that depends on the radii of 
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the stars relative to the orbital separation, H is the atmo- 
spheric depth at which currents cross magnetic field lines 
and return back to the non- magnetic NS, and Ad is the 
thickness of the arc-like cross section of the current layer 
in the magnetic NS's atmosphere. Given the uncertain- 
ties in these parameters, I approximate J ~ H/ Ad ~ 1 
for these initial estimates. This results in 



fti = 1.2 x IQ~ M K 



Hi 



10 2 B- 



\ (13) 



for the resistivity of the magnetic NS. 

The other key resistivity is that of the space between 
the NSs, 7^ sp acc . Given the large potential estimated in 
equation (flOl) . charged particles are likely accelerated to 
sufficient energies for pair production. This creation, an- 
nihilation, and the subsequent radiation, coupled with 
other dissipative processes like Alfven wave generation 
and damping, implies that this space may dominate 
the resistivity of the circuit (for example, see the dis- 
cussions of dissipation in m agnetar magnetospheres by 
iThompson fc Duncanl[l9 95) . In the relativistic limit, one 
could approximate 



4tt/c = 4.2 x 10 



10 s cm" 1 , 



(14) 



which is over ten orders of magnitude larger than the 
estimated resistivity of the magnetic NS ( e g . [T51 ) . This 
is the value used by iMc Williams fc Levin! (|2011 ) for the 
resistivity of a black hole event horizon (also see lDamourl 
H98l IThorne et al.lll986ft . But given the complications 
of what is potentially occurring in the magnetosphere, 
it is not at all clear that such an extreme resistivity is 
accurate. So in the interest of understanding the full 
range of potential solutions, I will explore 7?. sp aco as a 
free parameter. 

The dissipation can therefore take very different forms 
depending on the relative resistivities. When the mag- 
netic NS dominates (Hi 3> H spaC c), most of the dissi- 
pation occurs within the surface of the star where the 
currents run (see Figure [TJ), mainly producing heat with 
a rate of 



E 1 « 2.4 x 10 47 a4#16 



10 2 s- 



17/3 



erg s 



(15) 



This result shows that energy dissipation from magnetic 
effects is potentially large, but it depends strongly on how 
asynchronous the magnetic NS is (via the ratio ui/ui). 
This is turn is mediated by the competing effects of grav- 
itational wave emission which makes the binary more 
asynchronous and magnetic torques, which try to syn- 
chronize the magnetic NS. I explore this competition in 
more detail in the next section. 

A potentially interesting issue is that as the dissipation 
heats the surface of the magnetic NS, it would increase 
its conductivity, which would in turn increase the dissi- 
pation, and so on. A more detailed investigation of the 
surface layer structure is needed to assess what kind of 
feedback loop could potentially result. Given these un- 
certainties, I keep K fixed with orbital separation with 
the caveat that it could in principle change as the binary 
coalesces. 



mainly occurs in the space between the stars with a value 

4 P acc « 6.8 x 10 36 /4 (l^i) 14/3 B_1 - 

(16) 

A resistivity this high implies a fairly weak observational 
signal, at least for typical pulsar-strength fields. If the 
magnetic moment is increased to fj, rj 10 33 G cm 3 , as is 
relevant for magnetars, and the binary is close to merger, 
then this dissipation can be greater than ~ 10 45 erg s . 
This is consistent with the prefactor and seari ngs found 
in equation (5) of Hansen fc Lyutikovl (|2001( ). which I 
explore further in H4.'2{ So even in the limit of the 
strongest possible resistance, an observational signature 
of a merger is possible for magnetar strength fields. 

Finally, there is an additional source of dissipation that 
is not associated with the circuit itself, but should be 
mentioned. Since the non-magnetic NS is effectively a 
perfect conductor, the magnetic field lines cannot pene- 
trate this star and must be reoriented around it, which 
requires work. Alternatively, one can also think of this 
in terms of the current that is generated to produce a 
magnetic dipole in the opposite orientation to the verti- 
cal magnetic field (as was discussed in §2.1(1 . The rate 
of work done to move the field lines is roughly just the 
flux of energy density in magnetic field lines that are 
intercepted by the orbit of the non-magnetic NS, 

,, 2 ,,2 p2, ,13/3 _ 

E mh « x nR\a{uj - Oi) - 



16(GM) 5 / 3 u ' 

(17) 

Although this is potentially comparable to the dissipa- 
tion in the space between the NSs (eq. [15]). it is less 
important because it does not have a corresponding ob- 
servational signal. It does sap energy from the orbit, 
but this is a negligible effect in comparison to the grav- 
itational wave emission. Since this will not appreciably 
alter the inspiral, I ignore it for the remainder of this 
work. 

3. MAGNETIC TORQUES AND DISSIPATION ESTIMATES 

As shown in Figure [TJ the interaction of the current 
across the non-magnetic NS with the dipole field provides 
a torque on the magnetic NS, Nb, which spins it up 
when uj > Q\. The total magnetic torque from both 
hemispheres is 



N B = 



4fiR 2 l 



4$^ 



cri (Hi + 2K : 



(18) 



space ) 



In the limit that Hi 3> H sp ncc, Nb ~ 2 ja-\ . S imilar 
expressions t o this are derived in IWu et al. I (|2002|) and 
iDong fc Lml (|2004h . Note that N B oc ai oc u - Cli, so 
that indeed Nb > (spinning up the magnetic NS) when 
u> > Cli. 

The orbital angular momentum of the binary, which is 
given by J orb = (Ga/M) l / 2 MiM 2 , evolves as 



In the limit when 1Z S 



Att/c^> IZi, the dissipation 
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Np 



(19) 
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where J gw is given in equation (Q}. I ignore any tidal 
effect s since they cannot sign ificantly synchronize the bi- 
nary (Bild sten fc Cutlerlll992f) . For solid-body rotation, 
the spin of the magnetic NS evolves as 



dt 



(/ifii) « N B + N p , 



where 



6c 3 



(20) 



(21) 



is the torque from pulsar dipole spindown. 

Before solving the time evolution of the system nu- 
merically in [21 it is helpful to consider some analytic 
estimates that highlight the principal effects introduced 
by the inclusion of magnetic torques. There are a few key 
timescales that will be important in determining the evo- 
lution of the system. The first is the gravitational wave 
timescale, which was already presented in equation (J5]> . 
Another key timescale is that associated with the mag- 
netic torquing 



tb = 



N B 

27/i 31 2 ^r 6 1 ( 



4$ 2 \uj) \ 

-11/3 



2TZ, 



space 



10 2 s" 



211. 



Hi 



(22) 



Finally there is the timescale for spindown from magnetic 
dipole emission, 



fliii 

\NJ 



-6.7 x 10 3 /xii 2 



10 2 s- 



yrs. 



(23) 



Even if the magnetic NS is nearly synchronized to the 
orbit, the dipole spindown is negligible, so I ignore it for 
the remainder of this study. 

The evolution of the system is subject to the compet- 
ing effects of the gravitational wave emission promot- 
ing asynchronicity (thus increasing a± ) and the magnetic 
torques trying to synchronize the magnetic NS (forcing 
(Ti toward zero). This drives the system toward an equi- 
librium where dai/dt « (as was found for the compe- 
tition between tidal tor ques and gravitational waves in 
white dwarf binaries by lPiroll20l"T ) . The time derivative 
of the tidal forcing frequency is 



dai 
~dt 



= 2 



did 

~dT 



dtti 
~dT 



= 2 



i gw 



TB 



(24) 



Setting dai/dt w 0, 1 estimate the steady-state tidal forc- 
ing frequency to be 



11. 



TB 

T gw 



2K S 



72-1 



.(25) 



For a non-spinning NS, the ratio would be ui/uj = 2, 
so this shows that the magnetic torques act to make the 
magnetic NS more synchronous with the orbit. 



This result only holds when tb < 2r gw , and therefore 
magnetic torques are only important if the field exceeds 
a critical strength of 



B crit = 1.8 x 10 12 X 



ie 1/2 ( 



10 2 s- 



-1/2 



1 



211, 



space 



Tlx 



1/2 



G. 



(26) 



This critical field becomes smaller with higher orbital fre- 
quencies (.Bcrit oc u;^ 1 / 2 ), showing that magnetic torques 
increasingly play a role as the NSs get closer to merger. 

Whether or not magnetic torques are important, there 
will be a large electrical power dissipation on the mag- 
netic NS as long as 72. spacc is negligible. At early times 
when B < J5 cr i t , I substitute <t\/uj « 2 into equation (TTT 
to find 



£a « 9.6 x 10 4 Vli^i6 (^-t) 



17/3 



erg s , 
B<B crit (27) 

When B > B cl -i t , substituting equation (|25|) into equa- 
tion (|15l) results in 

E x « 8.7 x lO^si^ie 1 (l^i) 1173 ^ 

B > B ciit (28) 

where I have assumed in both cases that 72. spaco IZi. 
The dissipation therefore is expected to follow a broken 
power law, with Eb oc w 17 / 3 at early times and Eb oc 
w 11 / 3 at late times. The reason for this flattening is 
that as the magnetic torque becomes more important at 
a closer orbital separation, the synchronization of the 
magnetic NS increases and <Ji/w decreases. This in turn 
decreases the induced e.m.f. 

4. NUMERICAL CALCULATIONS 

To explore the time evolution of the coalescing binary 
in more detail, I now turn to numerical integrations of 
equations (TIT))) and (j2"D)) using the prescriptions presented 
in the previous sections. These confirm many of my an- 
alytic estimates and provide more details about time- 
dependent effects. The integrations continue forward in 
time until th e stars reach the tidal disruption separation 
(|Kopa]|IT959L adding Fishbone's 1973 10 percent strong 
gravity correction) 



at 



2AR 2 



M 
W 2 



1/3 



(29) 



where I have assumed M 2 < M\. Unless noted otherwise, 
in all these calculations I use M\ = M 2 = 1.3M©, R\ — 
i?2 = 12 km, and fix K = 10 16 s _1 . This is because my 
main focus is understanding how varying the magnetic 
field and resistivity of the space between the NSs impact 
the evolution. 

4.1. Results when 7^ spa co <§C IZi 

I begin by exploring the case when 7?. S pacc <C Hi, so 
that the resistivity of the space between the stars is ef- 
fectively negligible. This makes a direction connection 
with the analytic results from the end of i|3] In the 
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Fig. 3. — Spin of the magnetic NS and the orbital period deriva- 
tive as a function of the gravitational wave frequency f = u)/ir 
when the resistivity of the space between the NSs is ignored 
(T^space -C IZi). The solid lines are labeled by the dipole mag- 
netic field strength (in units of Gauss). For all the calculations I 
use Mi = M 2 = 1.3M , R ± = R 2 = 12 km and fix K = 10 16 s" 1 . 
The dashed line marks the condition f2i = uj in which the magne- 
tized NS is synchronized with the orbit. Magnetic fields stronger 
than those plotted here simply push the magnetic NS even closer 
to synchronization. 



top panel of Figure [3l I plot the spin frequency of the 
magnetized NS as a function of the gravitational wave 
frequency / = oj/it. The solid lines are labeled by the 
dipole field strength, which varies from 3 x 10 10 G to 
10 13 G, from bottom to top. The dashed line is the rela- 
tion Oj = uj, showing that as the field increases the mag- 
netic torques drive the magnetic NS toward being syn- 
chronized. For a field of 10 12 G, synchronization starts 
becoming important around a frequency of / « 200 Hz. 
This is roughly consistent with B cr it estimated from plug- 
ging uj = 7r/ w 600 s" 1 into equation (ffij]). 

As the magnetized NS spins up, it could in princi- 
ple experience rotational instabilities which would result 
in a triaxial shape and produce additional gravitational 
waves. These instabilities scale with the rotation param- 
eter P = Erat/lW], where E TOt is given in equation (Ht 
and I use the prescription given in Lattimer fc Prakashl 
(pOOll) for 



\W\ w 0.6Mic 2 



GMx/Rxc 2 



1 - 0.5(GMi/i?ic 2 ) 



(30) 



Dynamical bar-mode instabilities occur for /3 > 0.27 
(|Chandrasekhan1l969l ). and secular instabilities for /3 > 
0.14, d riven by gravitation al radiation reaction or vis- 
cosity (|Lai fc Shapiro! ^995). Even for the latter case, it 
requires / > 2.8 x 10 3 Hz for instability, so in practice 
the binary always merges before sufficiently high spin 
frequencies can be reached. 
In the bottom panel of Figure [3l I plot the fractional 
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Fig. 4. — The difference in the number of orbits per logarithm 
orbital frequency due to magnetic torques when 7£ S pace -C TJi. 
Each solid line is labeled by the magnetic field strength (in units 
of Gauss) . All other quantities are fixed the same as in Figure [3] 
From this presentation, one can see at what frequencies the most 
number of orbits are being lost. For a magnetic field > 6 X 10 12 G, 
no additional orbits are lost because the magnetic NS is already 
effectively tidally locked. 



difference of the orbital period derivative 

2 / „,r \ 3 



\P\ = 6tt 



fJc 



t-b 



\GM J \MxM 2 



M 



orb I 



(31) 



in comparison to that given from just gravitational wave 
emission |P gw | (which is found by just substituting | J gw 

for |J or b| hi the previous expression). In general, be- 
cause magnetic torques remove angular momentum from 
the orbit and put it into the magnetized NS, the absolute 
value of the orbital period derivative increases (i.e., be- 
comes more negative). The maximum effect occurs when 
the magnetized NS is synchronized with the orbit. The 
bottom panel of Figure [3] shows that the period deriva- 
tive is larger by as much as « 10 — 40% at late times. 
It should be emphasized that this is the maximum pos- 
sible deviation. For magnetic fields > 10 13 G, the period 
derivative cannot change by any more because the mag- 
netic NS is already nearly synchronized with the orbit 
at large separation. The larger period derivative means 
that the binary experiences less orbits while its inspi- 
ral is at frequencies detectable by ground-based interfer- 
ometers. This will create difficulties in using matched- 
filtering techniques to pull these binary inspiral signals 
out of the noise and making re liable parameter estima- 
tions (ICutler fc FlanaganlH99l . 

Although the period derivative changes by the largest 
amount at high frequencies, there are many more orbits 
at low frequencies. This means that the largest differ- 
ence in the number of orbits (which is what is important 
for matched filtering) may be at earlier times. To bet- 
ter quantify how many fewer orbits are expected and at 
what times, I plot the differential number of orbits N per 
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Fig. 5. — The ratio of the tidal forcing frequency to orbital fre- 
quency, ai/cj, as a function of the gravitational wave frequency 
/ = uj/tt (solid lines). Each curve is labeled by the magnetic field 
strength (in units of Gauss), with all other quantities kept fixed 
with the same values as in Figure [3] The dashed curves show the 
ratio of the magnetic torquing timescale to the gravitational wave 
times, ra/rgw, for each corresponding magnetic field, which is pre- 
dicted to be a good estimate for the ratio 0"i/oj and be proportional 
to u)^ 1 in equation H25I I. 

logarithm frequency in Figure 21 where 
dN 



din/ 



= P 



(32) 



and similarly dN gw /dln f = \Pgw\~ 1 . Such a plot shows 
at what orbital frequencies the most number of orbits 
are being lost. From Figure SJ one can see that as the 
magnetic field increases, the majority of the orbits are 
lost at earlier times during the inspiral. For example, at 
a field strength of 10 12 G, a few orbits are lost around a 
frequency of / w 200 Hz. Such differences in the inspiral 
will have to be accounted for in a matched-filtering anal- 
ysis (which would of course also include post-Newtonian 
effects that I have omitted for simplicity). 

I next check the relation estimated in equation (|25[) 
for the tidal forcing frequency in Figure [5] In each case 
at early times o~i/lj ~ 2, which corresponds to a low 
spin frequency for the magnetized NS. But as the mag- 
netic torque increases at high orbital frequency, ui/uj 
roughly asymptotes to the ratio of the magnetic torquing 
timescale to the gravitational wave timescale, Ts/r gw , 
which is plotted as a dashed line. At magnetic field 
strength < 10 12 G the magnetic torques are less impor- 
tant and equation (|25|) no longer applies. 

In Figure I calculate the electrical power dissipation 
on the surface of the magnetic NS, Ei, for a range of 
magnetic field strengths. In the top panel I plot this as 
a function of /. This shows a dependence on the fre- 
quency that roughly matches the scalings given in equa- 
tions (|2"T)) and (|2"51) at low and high fields, respectively. 
Due to the effect of magnetic torques, the highly mag- 
netic NSs show more dissipation at low frequencies and 
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Fig. 6. — The electrical dissipation power plotted as a function 
of both the gravitational wave frequency / = ui/tt (top panel) and 
the time until merger (a at, bottom panel). The different line 
types correspond to different magnetic fields, as indicated by the 
key in the upper panel. 



low magnetic field NSs show more dissipation at high 
frequencies. Although the power can get rather large at 
Eb ~ 10 49 — 10 52 erg s _1 , this only last for a short time. 
In the bottom panel of Figure [HI I plot these same en- 
ergy dissipation rates as a function of time until merger 
(a w at). This shows that a rate of > 10 49 erg s _1 is 
only possible for the last ~ 1 s. Such energetics and 
timescales are similar to short GRBs, a fact I discuss in 
more detail in the conclusion. 



4.2. Results when TZ S 



is Non-negligible 



The largest uncertainty in these calculations is the re- 
sistivity in the space between the stars lZ space . Thus far 
I have been assuming that 7^ S pacc is negligible, but in- 
teresting phenomena potentially occur no matter what 
value it takes. In Figure I fix the magnetic field at 
10 13 G and explore the effect of increasing 7^-spacc on the 
dissipation from the magnetic NS E\ (top panel) , and the 
dissipation in the space between the stars -E S pace (bottom 
panel) . Each line in Figure [7J corresponding to a differ- 
ent value of 7^-spaco in units of the maximum resistivity 
47r/c according to the key in the bottom panel. This 
shows that there is some give and take between these 
two rates of dissipation, and independent of the exact 
value, it appears that the dissipation somewhere should 
be observationally important. 

On the other hand, the deviation of the coalescence 
period derivative from just a point-particle inspiral is 
strongly effected by 7?. S pacc, as is show in Figure [3 
The top panel shows the fractional change of the pe- 
riod derivative. When 7Z spacc is small, I effectively get 
the same large deviation of the period derivative as was 
found in Figure [3] before. These differences persists as 
long as 7?. S p aco < IZi. But once these are about equal 
when lZ spaC c < 10~ 10 x 47r/c, the period derivative is ex- 
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Fig. 7. — The dissipation rates on the surface of the magnetic NS 
(El, top panel) and in the space between the NSs (-E spa cc, bottom 
panel) as 7?. spa cc changes. The value of 7£. S pace are in units of 47r/c 
as denoted by the key in the bottom panel. For all the calculations 
I fix B = 10 13 G, and the masses and radii are the same as in 
Figure [3] 
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Fig. 8. — The period derivative (top panel) and number of orbits 
lost (bottom panel), for different values of 7?. S pacc in units of 4ir/c, 
as labeled. This demonstrates that the differences in the period 
derivative and number of orbits are strongly dependent on "R. S pacc- 
For all the calculations I fix B = 10 13 G, and the masses and radii 
are the same as in Figure [3] 

pectcd to be consistent with point-particle inspiral. In 
the bottom panel of Figure [5J I better quantify this in 
terms of the number of orbits that are lost at a given 
frequency (similar to what was done before in Figure [4]). 
In all cases the peak dissipation rates are reached near 




Fig. 9. — The dissipation rates at the moment T g 



0.1 s which 



corresponds to / m 480 Hz, from both the surface of the magnetic 
NS [Ei, top panel) and in the space between the NSs (-Espace, bot- 
tom panel) as 'R.apace changes. For all the calculations the masses 
and radii are the same as in Figure [3] 

the end of the NS inspiral, and this provides the most 
exciting opportunity for a transient signature coincident 
with gravitational wave emission. To better quantify 
what is expect, I plot the dissipation rates at the mo- 
ment r gw = 0.1 s in Figure [S] This corresponds to an 
orbital frequency of oj = 1.5 X 10 3 s _1 or a gravitational 
wave frequency of / rj 480 Hz. This figure again displays 
the give and take between the two dissipation locations 
as 7£spacc changes. It also highlights that the most po- 



tentially optimistic case comes when 7Z S 



Tl x . Also 



note that in the limit 7^. spaco = Att/c and B = 10 15 G, 

the dissipation in the spa ce between the NSs fipace is 
similar to the estimates in [Hansen fc Lvutikovl (|2001[) . 

5. DISCUSSION AND CONCLUSIONS 

I investigated the inspiral of NS binaries when one NS 
has a magnetic field. Using a unipolar inductor model, I 
assessed what typed of currents and associated electrical 
dissipation rates are expected. The main uncertainty 
in this analysis was the resistivity between the stars, 
T^space, but for whatever value it takes in nature there 
are interesting consequences to explore. When 7?. S pacc 
is large, electrical dissipation comes predominantly from 
the space between the NSs. This implies specific obser- 
vational predictions, which I discuss below. When 7?. S pacc 
is small, the competing effects of magnetic torques spin- 
ning up the magnetized NS and gravitational wave emis- 
sion causing the binary to coalesce result in a steady- 
state asynchronicity. I explored this both analytically 
and with numerical integrations of the binary time evo- 
lution. The loss of angular momentum from the binary 
that is used to spin up the magnetic NS impacts the dy- 
namics of the inspiral. There is also a large associated 
dissipation of currents in the surface of the magnetic NS 
that may again have observational consequences. 
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5.1. Observational Signatures 

When 7?. S p acc is non-negligible, the dissipation mainly 
occurs in the space between the stars. This will poten- 
tially result in a couple o f different signatures , which are 
analogous to the effects Hanse n~fc Lvutikovl (|2001[ ) de- 
scribe in some detail, but I quickly summarize the main 
results here for completeness. The electric field gener- 
ated around the non-magnetic NS accelerates particles to 
sufficient energies for pair production, analogous to radio 
pulsars. In pulsars, the energy extraction is limited to the 
polar caps, but for the non-magnetic NS, all field lines ex- 
tend away, so that the po lar cap effectively en c ompa sses 
the whole star. From this lHansen fc Lvutikovl (|2001h es- 
timate that for a magnetar a radio flux of ~ few mJy 
is possible. On the other hand, there are many uncer- 
tainties in whether the rad io emission will suppre s sed by 
the s trong magnetic field ()Usov fc Melrose! [l996t lAronsI 
1998) a nd the shroud of p lasma within the magneto- 
sphere (|Goldreich fe Julian! fl969l) . Also, it is uncertain 
what kind of delay interstellar d ispersion may create 
(|PaImerl ll993: Lipu nova et al.lll997l ). causing the radio to 
follow the actual merger by hours or more. Potentially 
more promising in the case of a large lZ sp&cc is the energy 
dissipation of accelerated particles and Alfven waves in 
the space between the stars. This would create a hard 
X-ray precursor from a relativistically expanding wind 
of pairs and photons that would occur ~ seconds before 
merger. 

In the limit when 7^ spa cc is small, the majority of the 
dissipation occurs in currents in the surface of the mag- 
netic NS. This can reach rates of > 10 49 erg s _1 during 
the final ~ 1 s before merger, as shown in Figure |H1 Al- 
though the energetics and timescale of the dissipation are 
similar to short GRBs, it is not clear what the observa- 
tional outcome is here. The majority of the dissipation 
will take place where the field lines connecting to the 
two NSs intersect the magnetized NS. This presumably 
is released as heat. If a small fraction of mass is ablated 
from this energy injection, it could produce a pair fire- 
ball that may be like a GRB. The main difficulty with 
making this connection is that it is not clear what col- 
limates the flow sufficiently to match the observed jets 
(although there is some evidence that at least some short 
GRBs have relatively wide jets, IGrupe et al.M2006h . but 
this clearly merits a more detailed study. 

5.2. Future Work 

The effect of magnetic torques on binary NS inspi- 
ral dynamics deserves further study, especially since 
matched-filtering is anticipated to play an important role 
in detecting these events with ground-based interferom- 
eters. When the inspiral is altered by magnetic torques, 
it will likely lead to misestimations of binary parame- 
ters. Future work should quantify which parameters are 
most impacted and what are the sizes of the induced 
uncertainties. In addition, compact binary coalescences 
have historically been identified as promising sites for 
studying the strong-field dynamics of general relativity, 
and potentiall y testing for other theories of gravity (see 
iLi et al.l 120121 and references therein). Magnetic fields 
will make such efforts more difficult. 

On the other hand, given some of the electromag- 
netic predictions from these magnetic interactions, there 



may be correlations between different inspiral signals 
and X-ray precursors (~ sec before merger) and/or 
later radio emission (~ hrs after). This adds to a 
host of other potential electromagnetic phenomena that 
have been theorized to be associated with compacter 
mergers , including crust c racking in the seconds before 
merger (iTsang et al.ll2012f), kilonovae in the day follow- 
ing merger ( )Metzger et al.l 120081 I2010D, and r adio tran- 
sients a few weeks later ()Nakar fe Piranl 120111) . Each of 
these emission mechanisms probe different, complemen- 
tary aspects of the merging binary, and when taken to- 
gether could provide a detailed picture of exactly what 
type of objects are being detected with gravitational 
waves. 

A potential complication is that throughout this study 
I have assumed that the less magnetic NS has a com- 
pletely negligible field. This is a good approximation at 
late times and for very strong magnetic fields for the mag- 
netic NS. But because a dipole field drops off so rapidly, it 
does not take a very large separation for each NS to have 
its own magnetosphere. For example, in the case of the 
binary pulsar PSR J0737-3039, only for a < QR X will the 
higher magnetic field really dominate. Prior to this time, 
a different analysis mi ght be required that takes into ac- 
count both fields (e.g . iDemorest et al.ll2004l : iKaspi et al.l 
[2^0l:ILvutiko^l200llArons et alJl2005l^ 

Finally, another important problem is the coupling 
between the magnetic field, the magnetic NS's crust, 
and its core. In the present work I assume that the 
magnetic torque simply acts over the entire magnetic 
NS, and that it remains rigidly rotating. But since 
the field is primarily locked into the crust and the 
torques are concentrated at the magnetic footpoints (as 
shown in Figure [TJ, it is not clear if my assumptions 
hold in detail. For example, a solid NS crust has an 
associated shear modulus n cr ~ 10 30 erg cm -3 (substi- 
tuting typical values for the density and compo sition 
into the calculations by iStrohmaver et all [1991). In 
comparison, typical torques that we consider have 
values of N B « 2E B /a! ~ 10 46 - 10 50 erg. Even 
when averaged over the entire NS crust, with a typical 
thickness of ~ 0.1R, this gives a torque density of 
~ 10 29 — 10 33 erg cm -3 . Given that the torque is 
actually exerted on smaller regions of the crust, it seems 
likely that they may exc eed the crust breaking strain 
([Horowitz fc Kadaul ' l2009T ). A more detailed study of 
the crustal coupling is therefore needed. This will help 
develop a better understanding of how the torques are 
transmitted through the NS. If crust cracking occurs, 
it may lead to interesting precurso r phenomena, suc h 
as in the observations reported by (jTroia et alJ 12010). 
and suggested to be due to crustal breaking (albeit by a 
different mechanism) bv lTsang et al.l (|2012f) . 
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